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We have analyzed different correlation functions in a realistic spin-orbital model for half-doped 
manganites. Using a finite-temperature diagonalization technique the CE phase was found in the 
charge-ordered phase in the case of small antiferromagnetic interactions between t2 g electrons. It is 
shown that a key ingredient responsible for stabilization of the CE-type spin and orbital-ordered state 
is the cooperative Jahn- Teller (JT) interaction between next-nearest Mn +3 neighbors mediated by 
the breathing mode distortion of Mn +4 octahedra and displacements of Mn +4 ions. The topological 
phase factor in the Mn-Mn hopping leading to gap formation in one- dimensional models for the CE 
phase as well as the nearest neighbor JT coupling are not able to produce the zigzag chains typical 
for the CE phase in our model. 

PACS numbers: 75.30.Kz, 71.70.Ej, 75.30.Et, 75.10.-b 



I. INTRODUCTION 

Half-doped perovskite manganites Rex/2 Ai /2M11O3 
(Re-rare earth, A-alkaline earth) exhibit very specific 
properties quite different from manganites with other e g 
electron concentrations i& Generally, for large size Re 
and A cations or when external pressure is applied 3 these 
compounds are in the A-type antiferromagnetic (AF), 
metallic state at low temperatures showing no or only 



experiments are interpreted in terms of almost perfect 
CQ 6 ^ while others^ 3 , are considered to be consistent 
with a smooth charge density wave (CDW) of the e g 
electron density. Although, it is plausible that the 
complex spin and orbital arrangement in the charge- 
ordered manganites can result from the competition 
between DE24 and AF superexchange interactions (~ 
Jap), the role played by the CO and Jahn- Teller (JT) 
interactions22iSi*22iS in the system is still unclear. Fi- 



weak sign of charge ordering, and the occupied e g elec- nall y * system* l^e e.g. Pr 1/2 Ca 1/2 Mn0 3 CO and CE- 



tron states are predominantly x — y like. Recently it has 
been reported that the A-phase may develop an intrinsic 
charge-stripe modulation, which controls the transport 
properties^ With decreasing size of the cations, which 
implies a decrease of the bandwidth, charge ordering is 
observed which at lower temperatures is accompanied 
by the formation of peculiar ferromagnetic (FM) zigzag 
chains 5 - 6 ' 7 that are staggered antiferromagnetically. This 
CE spin order, a notation introduced in the pioneering 
work of Wollan and Koehler— and Goodenough-;, is ac- 
companied by a checkerboard charge order (CO) and di- 
rected occupied orbitals on the Mn 3+ sites amplifying 
double exchange (DE) along one lattice diagonal. The 
CE-structure is a particularly fascinating manifestation 
of the control of magnetic order due to the orbital de- 
gree of freedomiSiii in combination with charge order- 
ing. Although a number of approaches have been fol- 
lowed to explain this structure, the mechanism for its 
stability is still not clear and hotly debated with em- 
phasis put on the role of the DE t 12 i 13 i 14 the Coulomb 
interactionji 5 -i§ii 7 ^i£ or the coupling to the lattice de- 
grees of freedom, 19 i 20 i 21 i 22 i 23 i 24 i 25 

The CE phase has been observed experimentally both 
in cubic (Nd,Pr) 1/2 (Sr,Ca) 1/2 MnOj 2 i 7 i 26 i 27 i 28 and in lay- 
ered LaSr 2 Mn 2 0-^ and Lai/ 2 Sr 3 / 2 Mn04^* 6 " manganites. 
In some cases it can coexist with the A-type AF spin 
ordering. 29,30,31 Regarding CO in the CE phase, some 



type orbital correlations are found to develop well above 
the Neel temperature 3 ^ 7 , or even coexist with FM spin 
stated which indicates that these orbital correlations are 
not primarily driven by magnetic interactions. 

This work focuses on a microscopic understanding of 
the stability of the CE-phase which appears to be the 
ground state of narrow band manganites at and near 
quarter-filling. We argue here that the key interaction, 
apart from the nearest-neighbor Coulomb repulsion V 
which supports the charge ordering^ is the next-nearest 
neighbor JT interaction. This interaction between or- 
bitals at next-nearest neighbor Mn 3+ ions emerges from 
the JT-distortions of the Mn 3+ octahedra which is ampli- 
fied by the breathing distortion of the intermediate Mn 4+ 
octahedra. We show that this interaction generates a 
narrow regime in the V-Jaf phase diagram where the 
CE-phase (and sometimes the C-phase) is stable. Other- 
wise one encounters a homogeneous FM phase at small 
Jaf fully controlled by DE and a conventional nearest- 
neighbor AF phase at larger Jaf- 

An important feature of further neighbor JT interac- 
tions are the displacements of the Mn 4+ ions. We are 
only aware of work by Radaelli et a/4°» where detailed lat- 
tice coordinates for the CE-phase in Lai/ 2 Cai/ 2 Mn03 are 
reported. In particular they observed that the 4 planar 
0(2) move towards the Mn 4+ ion and the Mn 4+ ions are 
found to be displaced from their regular perovskite posi- 
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tions, while the Mn 3+ ions are not shifted. This observa- 
tion confirms the importance of the breathing distortion 
and also shows the mechanism of the further neighbor 
JT-distortion in the CE-phase. As a consequence the ap- 
pearance of new Bragg reflections in the CE-phase would 
be associated not primarily with quadrupolar electronic 
orbital order but rather with the displacements of ions 
associated with the CE orbital order 4£ 

An alternative long-ranged JT based interaction was 
recently proposed by Khomskii and Kugel^i and by 
Calderon, Millis and Ahn4£ They employ elasticity the- 
ory to derive the orbital interaction between JT-distorted 
Mn 3+ C>6 octahedra for the perfectly charge ordered case. 
We have also investigated the consequences of this inter- 
action. Our numerical simulations confirm the estimate 
presented in Ref. ^] that for perfect charge ordering the 
AF CE phase is stable. Yet in this model this requires 
large nearest neighbor Coulomb repulsion V, while for 
moderate values of V the CE phase gets unstable. 

The finite temperature diagonalizationi&iM4 used in 
our study allows to include the full interplay of spin 
and orbital states, and to monitor the onset of differ- 
ent orders as function of temperature. This gives us the 
ability for an unbiased investigation of the formation of 
different spin-orbital orderings emerging from all multi- 
electron configurations in the cluster. In this paper we 
concentrate on the stability and temperature dependence 
of spin, orbital, and charge ordered states in the two- 
dimensional (2D) realistic model for half-doped mangan- 
ites. In our microscopic spin-orbital model CO devel- 
ops as a result of nearest-neighbor Coulomb repulsion V 
and as a consequence of further neighbor JT interactions, 
while the spin-orbital order results from the competition 
between DE (i.e. kinetic energy), AF superexchange and 
further neighbor JT interactions. As shown in Sec. IIII Dl 
the latter interaction can lead to a CO state at low tem- 
peratures even in the absence of inter-site Coulomb re- 
pulsion. Furthermore we argue that a topological sign in 
the hopping amplitudes of e g orbitals 13 invoked in one- 
dimensional (ID) models for the CE phase is not suffi- 
cient to explain the formation of CE spin-orbital struc- 
ture in the 2D model studied by us. 

Further motivations to put effort in a better under- 
standing of the interplay of spin, charge, orbitals and 
lattice in the formation of the CE-phase are the follow- 
ing: It has been argued that the intrinsic mechanism that 
leads to colossal magnetoresistance are tendencies to- 
wards CE-type charge and orbital ordering in the vicinity 
of x = 0.54^ Furthermore a CE-insulator to FM-metal 
transition has been observed in relatively small magnetic 
field, which is accompanied by a major change of the 
optical conductivity on a surprisingly large energy scale 
(~ 1 eV)4£i2£ Such dramatic changes in relatively weak 
magnetic fields are a clear manifestation of the subtle in- 
terplay between spin, orbital and CO. The present study 
provides a basis for further investigations of the tempera- 
ture dependence of e.g. the optical conductivity and the 
study of spin and orbital excitations in the CE-phase, as 



well as the effect of doping into the quarter filled state. 

The paper is organized as follows. In Scc.[n]we present 
the model Hamiltonian, and describe in Sec. Ill Bl the 
finite-temperature Lanczos method. Sec. lIIII contains nu- 
merical results for different two-site correlation functions 
characterizing the charge, orbital and spin structure. 
These correlation functions are evaluated as functions of 
temperature and collected in a form of semi-quantitative 
phase diagrams in three different regimes: (i) including 
nearest-neighbor effective orbital-orbital (00) interac- 
tions (Sec. IIII B|) : (ii) assuming the elastic form of 00 
coupling at different distances (Sec. IIII C() : (iii) consid- 
ering the local form of next- and nearest-neighbor 00 
coupling (Sec. IIII D|) . while in Sec. IIII Fl the effect of 
charge stacking is discussed. The paper is summarized in 
Sec. IIVI In the Appendix the stability of the CE versus 
the C phase is analysed for a ID band model with lifted 
orbital degeneracy. 

II. MODEL HAMILTONIAN AND NUMERICAL 
PROCEDURE 

A. Effective spin-orbital model 

We consider the generic FM Kondo lattice model for 
manganites^^ with the Mn e g electrons coupled to the 
core tig spins via the Hund coupling Jh- The model is 
augmented by inter-site Coulomb repulsion V and coop- 
erative JT interactions and treated in the limit of infinite 
on-site Coulomb repulsion between two e g electrons on 
the same site (U — > oo): 

H = Hb anc i + H Kondo + Haf + Hy + HoO- (1) 

The first term describes the motion of the e g electrons, 

H band = - (4 c 4A-+ H - c -)' ( 2 ) 

with a constraint that allows only for Mn 4+ and Mn 3+ 
configurations, d^ CT = d^^l-rii^) Y[ <r '( 1 ~ n ^<x'): where 
all other alternative e g states are projected out. Here and 
in the following the index £ (a) denotes the e g orbital 
(spin) orthogonal to the £ (a) one, respectively. The 
hopping matrix elements depend on the basis chosen. For 
the basis {\x), \z)} with \x) ~ x 2 — y 2 and \z) ~ 3z 2 — r 2 
they are given by, 

}- t -( 3 tV *\ (3) 

where, +(— ) refer to the a(b) direction, respectively. The 
second term of Hamiltonian stands for the Hund in- 
teraction between e g electrons spin and the S — 3/2 ti g 
core spin Si, 

HKondo = — Jh ^ Si • d^a^o-'di^i, (4) 
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resulting in the parallel alignment of spins between ti g 
and e g electrons at each site in the limit Jh t which 
corresponds to the realistic situation in manganites. The 
Hband + HfCondo part of the total Hamiltonian alone rep- 
resents the DE mechanism and yields a FM ground state 
with | a;) orbitals being occupied which is favored by the 
kinetic energy in 2D model. Such a fully spin polarized 
case at finite U lea ds to the orbital t-J model as consid- 
ered in Refs. I43l44i Here, we are interested in the more 
complex case where both spin and orbital degrees of free- 
dom play an active role. Therefore, an AF (Jaf > 0) 
coupling between nearest-neighbor ti g spins is also in- 
corporated, 

-ffAF = Jaf E Si • Sj . (5) 

<ij> 

In our numerical studies we assume S = 1/2 for the core 
spins which decreases the resulting Hilbert space consid- 
erably. We have tested that this assumption gives similar 
results as in the 5 = 3/2 case when Jh is renormalized 48 
by the factor 3. Moreover, to promote the charge order- 
ing as observed experimentally in half-doped manganites 
we include the Coulomb repulsion (V > 0) between e g 
electrons on neighboring sites ^ 

H V = V E ^i^njfCT', (6) 

with rii^o- = dl^di^cr. This term reduces the probability 
of occupancy of neighboring sites by e g electrons. 

The form of the inter-site 00 interactions included in 
the -ffoo term of the Hamiltonian will be specified in the 
following section. 



the trace of the thermodynamic expectation value is per- 
formed by a Monte Carlo sampling leading to the follow- 
ing approximate expression: 

L M 

< A *> * if EE ex p(-^)(^i)(^i^io 5 (8) 

1=1 j=l 

where the first sum is performed over a restricted number 
L of random states |Z) being initial states in the Lanczos 
algorithm to generate a set {'ipA of eigenfunctions of H 
with respective eigenvalues {Ej}. M is the number of 
Lanczos functions in the expansion of a given state |Z). In 
the same manner the partition function is approximated 
by, 

L M 

Z^^exp(-^)|(!|^| 2 . (9) 

1=1 3=1 

It has been showni^ that the results get very accurate 
already for L, M <C N st . Even though we exploit the 
translational symmetry the calculations are restricted to 
small clusters (N — v° x v° is used in this paper) as 
only S* ot = S? subspaces can be treated separately, 
while a similar symmetry in the orbital sector does not 
exist. Although a V8 x V8 cluster is large enough to cap- 
ture the CE-type correlations, boundary conditions can 
influence further neighbor interactions. To assess such 
finite size effects quantitatively, a full diagonalization on 
larger clusters would be necessary. In the case of the CO 
CE phase this would involve at least a 4 x 4 site cluster, 
which is far beyond numerical realization. 

III. NUMERICAL RESULTS 



B. Finite-temperature Lanczos method 

To investigate the possible spin and orbital orderings in 
our spin-orbital model Q we have calculated the temper- 
ature dependence of different two-site correlation func- 
tions given by the operator expectation value: 

i N °t 

(An) = -: £>l exp(-/3ff)A R |n), (7) 

71=1 

of two-site operators, = l/NJ2i -Bi-Bi+R with B\ = 

m, Tf, s?, or S?m with m = Y,^a n m<?'^ Here > P = 

1/ksT, N st is the number of basis states {\n}} that 
span H, N number of sites in the cluster and Z = 
J2n=i( n \ ex P( — PH)\n) is the partition function. 

The calculation of (7J would, however, require a com- 
plete diagonalization of a N st x N st matrix of the Hamil- 
tonian. Therefore, we use for the calculation of (An) 
a generalization of the exact diagonalization technique 
developed by Jaklic and Prelovsek4S In this approach 



A. Parameters and magnetic phases 

In the numerical studies we assume S = 1/2 for 
the core spins and the limit of strong Hund's coupling 
Jh jt = 15. The numerical results are almost identical for 
all Jh ^ 5i in the whole range of temperatures considered 
by us. The Coulomb repulsion V/t and the AF-exchange 
Jaf A are f ree parameters, and the data is summerized 
in form of a phase diagram in the V — Jaf plane. The 
hopping parameter t is chosen as energy unit. We focus 
in this study on CO states and thus assume finite pos- 
itive values for V in all the calculations. In perovskite 
manganites the effective hopping t strongly depends on 
the lattice parameters and distortions e.g. increases with 
contraction of the Mn— O bond length and also depends 
on the Mn— O— Mn bond angled Thus, its value can be 
only roughly estimated as t = 0.2 - 0.7eVi 43 i 50 i 51 i 52 

A major difficulty when searching for spin structures 
like those indicated in Fig. ^ an d the same applies to 
charge and orbital order, is that in a cluster with periodic 
boundary conditions translational (and rotational) sym- 
metries are not broken. Therefore we have to investigate 
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FIG. 1: Schematic view of different spin patterns found in 
the 2D model in presence of checkerboard CO (half-doped 
case). Thick (thin) arrows represent states with total spin 
S = 2 (S = 3/2) with (without) an e g electron present at 
a given site, respectively. The shaded lines in C and CE 
phases indicate the direction of DE carrier propagation. The 
preferred directional e g orbitals occupied at sites S — 2 are 
parallel to the shaded paths. 



correlation functions and compare their features with the 
different spin configurations (presented in Fig. [TJ. 

We shall consider three different forms for the orbital 
interaction term -ffoo an d show that one has to go be- 
yond the nearest-neighbor interactions to obtain the CE- 
type correlations. We begin the description of our calcu- 
lations with a simplified model (without Hoo term) but 
including the crystal- field splitting of e g orbitals ex E z , 



H z = -E : 



(10) 



where T ; z 



n- lza ) is the z component of 



the pseudospin operator. This orbital splitting accounts 
for the elongation of octahedra along the c-axis in lay- 
ered compounds, which favors \z) over |x) orbital occu- 
pation, and also counterbalances the trend towards oc- 
cupation of \x) orbitals, which is strongly favored by the 
kinetic energy in the 2D-model. Although, the model 
without Hoo interactions can lead to the stability of the 
CE-type orbital pattern^ the magnetic structure was 
FM. In a wide range of parameters (E z /t < 0, V/t > 0, 
Jaf/£ > 0) no stabilization of the AF CE spin phase was 
found, indicating that the inter-chain coupling is quite 
strong and cannot be neglected in realistic models for 
half-doped CO manganites. Such a conclusion agrees this 
recent density-functional calculations showing a consid- 
erable band dispersion normal to the chains^ Thus, we 
argue that the topological effect considered in ID model 
approximations, which gave a simple explanation of the 
existence of the CE spin-orbital phase at half-doping, is 
not likely to be the decisive mechanism in a more realistic 
2D case investigated here. Our simulations with purely 
electronic interactions do not provide evidence for ID FM 
chains that are coupled antiferromagnetically to neighbor 




FIG. 2: Phase diagram indicating the character of the dom- 
inant spin correlations found at small temperatures T < zj . 
The calculations include the nearest-neighbor OO effective in- 
teractions (k — 0.2t). '*' indicates representative points for 
which the correlation functions are presented below. 



chains, which is the basic assumption that could justify 
such ID models. Therefore, in the following sections we 
concentrate on the electron-lattice couplingS2*2ii22i£iii4 
and present a mechanism stabilizing the CE spin/orbital 
state in manganites. 



B. The role of nearest-neighbor orbital interactions 

The first aim of our study of the spin-orbital model 
||TJ in half-doped manganites is to investigate the role of 
nearest-neighbor OO coupling played in CO state. There 
are two different mechanisms contributing to this inter- 
action: (i) the cooperative JT effect^ and (ii) superex- 
change interactions!^^ Neglecting more complex spin- 
orbital terms£i*££ both effects can be described by the fol- 
lowing simple term derived before for undoped LaMn03 
compound^ 



Hoo = 2k ^Ty, 

<ij) 



(11) 



where the two-site orbital operator, Ty, between nearest- 
neighbor Mn sites in the {\x}, \z)} basis, 

Tu = T?-T? + 3If Tf T %/3(ITT/ + TfT?), (12) 

is described in terms of pseudospin operators: Tl + = 
Y,<jd\ X!T di za , and XT = djj^dixa-- The prefactor of 
the mixed term cx \/3 is negative in the a direction and 
positive in the b direction. The nearest-neighbor OO 
term can lead to attraction of electrons on neighboring 
sites with different orbital orientations but having the 
Coulomb repulsion JSJ) included in the model we avoid 
phase separation. 

As shown in Fig.|21for large values of the superexchange 
interaction between core spins (zJaf ~ t) an AF spin 
ground state is realized while with decreasing Jaf the 
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FIG. 3: Temperature dependence of the charge (a), orbital 
(b), and spin (c) two-site correlation functions calculated for 
neighbors at different distances R for the phase with FM spin 
correlations. Parameters: V = t, k = 0.2t, Jaf = 0.051 
Jh = 151 



kinetic energy of the e g electrons (~ t) starts to play 
an active role via the DE mechanism^ breaking some of 
the AF bonds which leads to ID FM chains which form 
straight lines (C phase). Further decrease of Jaf yields 
FM correlations in both a and b directions. As the kinetic 
energy is controlled by the inter-site Coulomb repulsion 
V and suppressed in the limit of V S> t we find the AF 
region increasing in size with increasing V/t. 

We discuss now in more detail the temperature de- 
pendence of the different correlation functions. Start- 
ing from the FM region (sec Fig. EJ) we find the evo- 
lution of weak CO for T/t < 1 with nearest-neighbor 
(1,0) charge correlations smaller than the (1,1) and (2, 0) 
next-neighbor ones which is consistent with the alternat- 
ing order of Mn +3 and Mn +4 ions. The orbital corre- 
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FIG. 4: Temperature dependence of the correlations as in 
Fig. [jj] but for the phase with AF spin correlations. Parame- 
ters: V = 2t, k = 0.21 Jaf = 0.31 J H = 15t. 



lations are positive between all neighbors indicating the 
preference of x 2 — y 2 orbital occupancy due to the ki- 
netic energy which is quite large in the FM state. Al- 
though all spin correlation functions are positive, they 
remain highly anisotropic as T — > as a result of the 
strong competition between the DE mechanism and AF 
superexchange interaction between core spins. 

Next, we consider the AF spin state in the case of 
large Jaf and V (see Fig.0}. Increasing Coulomb repul- 
sion (V = 2t) we find more distinct and almost isotropic 
alternation of the e g charge which now sets in already 
at T/t w 2. Here, the orbital correlations are correlated 
with the onset of the AF spin order at T/t ~ 0.5. For 
T/t — » the spin-spin correlations are strong and repre- 
sent an almost isotropic AF spin state. 

For intermediate values of Jaf and moderate Coulomb 
repulsion we find the phase with predominantly C spin 



6 




0.04 



0.03 



A 0.02 



0.01 - 



-0.01 



(b) 




0.2 



0.4 



0.6 




FIG. 5: Temperature dependence of the charge (a), orbital 
(b), and spin-charge (c) correlation functions calculated for 
the phase with C spin correlations. Parameters: V = t, k = 
0.2t, Jaf = 0.2t, J H = 15*. 



correlations at T/t — > (see Fig. EJ). As presented in 
Fig. El for Jaf A — 0.2 an d V/t = 1 the charge corre- 
lations are similar to those obtained for the FM state 
shown in Fig. El (a) while the spin ordering changes dra- 
matically. Moreover for T/t — ► the orbital correlations 
are weaker than those found in the FM state discussed 
above indicating the orbital tendency towards in-plane 
directional- type order. 

A convenient way to distinguish between phases with 
C or CE and other spin correlations is to evaluate a com- 
bined spin-charge correlation function being restricted to 
sites occupied by e g electrons, (S( «iS , i z +R rii + R.). Ana- 
lyzing the localized limit presented in Fig. ^ ( or its by 
7r/2 rotated version) one can easily see that, for the C 
phase one should find (S*thS? +r th+r) > for |R| = 2 
and (S?niS[ +R n i+R ) < for |R| = y/2 while for the 



CE phase only (S?mS? +n n i+R ) < for |R| = 2 with 
any other spin-charge correlations remaining small. In 
the case presented in Fig. El (S* niS? +R ni + n) functions 
are negative (positive) for nearest-neighbors along diago- 
nal (next-nearest-neighbors along a/b direction), respec- 
tively, in agreement with the arguments presented above. 
Furthermore, nearest-neighbor spin-charge correlations 
are small as a result of averaging over different cluster 
orientations (rotated by angles ±7r/2, ±tt). 

It is straightforward to see why the C-structure is fa- 
vored by the nearest-neighbor JT interaction. The C- 
phase can be viewed as a realization of the alternat- 
ing orbital structure of LaMnOa but with a modulated 
charge density. Hence this structure is compatible with 
the nearest-neighbor JT interaction, whereas in the CE- 
structure this interaction is frustrated. 



C. Elastic interactions 

The role of local lattice distortions due to the JT ef- 
fect, which creates an anisotropic strain field decaying as 
~ i?~ 3 and thereby promotes long-range orbital order- 
ing, was considered by Khomskii (see Ref . l4lf) . Following 
these arguments, we treat the OO term as an elastic in- 
teraction which is active also between e g electrons on 
next-nearest-neighbor sites along a or b direction, 

H e i = 2n Ty + 2n' ^ Ty , (13) 

<ij> 

where ((ij)) refers to next nearest-neighbors along a 
and b direction with k' = k/8 which follows from the 
~ j|)~ 3 decay of elastic interactions^ The existence 
of longer-range OO interactions is justified by its cooper- 
ative lattice-driven nature where an orbital flip on a sin- 
gle Mn site distorts not only the neighboring oxygens but 
also positions of other manganese ions around it. Here, 
the JT coupling preferring the tilting of the orbitals from 
the uniform orbital state leads to an alternating orbital 
order (see Fig. EJ) and in this way compensates for the 
kinetic energy gain in uniform \x) orbital state favored 
in the 2D model. Here we neglect non-central elastic in- 
teractions between neighboring Mn sites in (11) direction 
as their strength does not scale with k' (see Ref. Elf) and 
thus would involve introduction of an additional param- 
eter in the model. 

The CE phase is favored as compared to the C phase 
by the OO coupling between second neighbors (~ k!) 
along a/b direction. However, this interaction is only 
active (not screened) when strong charge ordering takes 
place (large V/t case) and here the CE phase can persist 
even for small values of Jaf as long as the e g kinetic en- 
ergy is sufficiently suppressed by a large V (see Fig. |7J| . 
Thus, one has to assume large inter-site Coulomb repul- 
sion (where almost perfect Mn +3 /Mn +4 charge ordering 6 
is realized) and small superexchange interaction [see an 
example for V/t = 3 and J A fA = °- 04 in Fi S- El- In 
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FIG. 6: Schematic configuration of the charge-orbital ordered 
CE phase. The sites A and B indicate two orbital Mn 3+ 
sublattices while empty circles stand for the Mn 4+ ions. The 
dashed-line box indicates a typical orbitals orientation around 
a given Mn 4+ ion. 
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FIG. 7: Phase diagram as in Fig. |21 including the nearest- 
and next-neighbor elastic interactions (k — 0.2t, k' = ft/8). 
'*' indicates a representative point in the CE region for which 
the correlation functions are presented here. 



the {|x), \z)} basis, which we use in most of the numeri- 
cal evaluations, the spin and orbital correlation functions 
with R = (2, 0) become large as T — > while all the oth- 
ers remain very small in the whole range of temperatures. 
Such small orbital correlations imply either that the re- 
spective orbitals are uncorrelated or that they order in 
such a way, that the occupied orbitals are linear com- 
binations with almost equal amplitude expressed in the 
basis used for the measurement of correlation functions. 43 
This difficulty is not present in SU(2) symmetric models, 
and appears here as consequence of the cubic symme- 
try in the orbital sector. Therefore, the same functions 
were evaluated also in the {(|x) ± \z))/y/2} basis as pre- 
sented in the inset of Fig. 0(b). In the new orbital basis 
the function (T^Tf + , 2 Q is more pronounced and sat- 
urates for T/t ~ 0.1 at {T?T7- +{2fl) ) ~ -0.12 close to 
its maximum value —1/8 indicating almost perfect CE 
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FIG. 8: Temperature dependence of the two-site correlations 
as in Fig.|K|but for the phase with CE spin and orbital corre- 
lations promoted by elastic interactions (k = 0.2t, k' = n/8). 
Other parameters: V — 3t, Jap = 0.04t, Jh = 15i. The in- 
set in (b) shows the orbital correlations in the rotated basis: 
{{\x)±\z))/V2}. 



orbital structure (see Fig. EJ). Shorter range (TfT? +R ) 
correlations should become small when averaged over dif- 
ferently oriented clusters (rotated by angles ±7r/2) [see 
Fig. n (d)]. From Fig. ^ (d) one can see that a simi- 
lar pattern as found for the two-site orbital correlations 
should be found for the (5' i z niS' i z +R ni+R) functions. In 
our cluster in ideal CE state only electrons at distance 
|R| = 2 can contribute to the charge-spin correlations 
(being negative in this case) which do not vanish when 
averaged over different cluster orientations. As shown in 
Fig- El (c) the (5f niS , f +R ni + R) correlations obtained here 
are consistent with the above characterization of the CE 
pattern with a strong negative signal dominating at dis- 
tance |R| = 2. The existence of charge-ordered insulating 
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FIG. 9: Temperature dependence of the average value of the 
T* orbital operator (a) and the inverse of the staggered charge 
susceptibility Xcltagg (b) obtained for different values of the 
inter-site repulsion V/t in the presence of elastic interactions 
(k = 0.2t, k' — k/8). Other parameters: Jaf = 0.041, Jh = 
15t. 



state concomitant with orbital and spin order of the CE- 
typeS** was observed in different half-doped manganites 
with small one-electron bandwidth. 1 

The average change in the orbital occupation can be 
studied calculating (Tf) which changes between +1/4 
and —1/4 for the uniform |a;) and \z) orbital states, 
respectively, for the case with one e g electron per two 
Mn sites. Furthermore, when the orbital state involves 
\x) and \z) orbitals in the same proportion one finds 
(T; 2 ) ~ 0. In Fig. [S] (a) we see that increasing the 
Coulomb interaction from V = t to V = 2t the orbital or- 
dering hardly changes for T/t — > with (T?) only slightly 
decreasing to (T?) ~ 0.13 for V = 2t. This indicates that 
as long as the system is in the FM state (see Fig. [7J the 
ground state is dominated by the x 2 — y 2 orbital pro- 
moted by the kinetic energy term. Further increase in 
the Coulomb inter-site repulsion to V — 'it leads to a 
drastic change in the orbital order when system enters 
the CE phase region (see Fig. 0). Here, (T?) ~ 0.03 
for T/t — ► indicating dominating role played by the 
(\x) ± \z))/y/2 orbitals. This fact is consistent with our 
results for the two-site orbital correlation functions [see 
inset in Fig. Mb)]. 

A change in the inter-site Coulomb interaction is di- 
rectly connected with the robustness of the charge or- 
dering as seen in Figs. 03- (a) and Fig. |S] (a). Al- 
though the characteristic temperature of the CO (Too) 
can be roughly estimated from the charge-charge corre- 
lation functions it can be more precisely extracted from 



1 

Ui -2 



c.stagg ■ 



3 2 



(14) 



which at high temperatures follows a Curie- Weiss law 
XZltagg tx (T — T C o) ■ From Fig. (b) one can easily 
estimate T co /t w 0.6, 1.4, and 2.4 for V/t = 1.0, 2.0, 
and 3.0, respectively, in agreement with the mean-field 
solution 56 " which for V ^> t predicts, Too ~ zV/4. 



D. Microscopic model for Mn — Mn JT coupling 

Although the model considered in the previous sec- 
tion with longer-range OO elastic coupling leads to the 
CE-type spin and orbital correlations, this phase is found 
only in a very extreme case of almost perfect Mn 4+ /Mn 3+ 
CO (that is for strong Coulomb repulsion) which may 
not be very realistic in the light of some experimental 
facts indicating the CDW character of charge redistri- 
bution between different manganese ions^* Therefore in 
this section we derive the further neighbor OO interac- 
tions from a microscopic model including apart from the 
JT-distortion of the Mn 3+ octahedra the breathing dis- 
tortion around the Mn 4+ ions as well. We then show that 
this interaction is capable to produce the zigzag (spin 
and orbital) phase in a region of small Coulomb repul- 
sion changing the previous phase diagrams in a drastic 
way. 

To assess the gain in lattice energy connected with the 
formation of the CE long-range orbital order we follow 
the work of Millis for LaMnOs by treating the lattice 
distortions classicallyi 5 ** 5 ^ The effect of lattice displace- 
ments is described by assuming arbitrary values of four 
different lattice distances di (i — 1,...,4) (see Fig. HOf) 
which are expressed by S x , 6' x - uniform deformation 
along the a and b direction of the Mn lattice, and Ui 
- displacements of O ions along the Mn-O-Mn bonds 
within the (a, b) planes: 



(15a) 
(15b) 
(15c) 
(15d) 



where b is the lattice constant of the ideal perovskiteiS 
Moreover, we assume e g charge stacking in the c direction 
with the respective Mn-Mn distance being 6(1 + S z ). 

In order to define orbital sublattices (A and B) we 
introduce a transformation which describes the tilting of 
orbitals by making two different transformations at both 
sublattices!* 5 - 
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FIG. 10: Schematic representation of the lattice distortions 
contributing to the second neighbor JT-interaction. The cor- 
responding CE orbital order for the 5 atom structure above 
is indicated in Fig.|S|by dashed lines. Different Mn 3+ , Mn 4+ , 
and O ions are presented by full, shaded and empty circles, 
respectively. The distances between the Mn 4+ ion and the 
Mn 3+ ions in (10) [(10)] and (01) [(01)] directions, respec- 
tively, are identical. The arrow indicates the displacement of 
the central Mn 4+ ion when d\ + d.2 7^ d L i + d,4. 



> 





AF 


* CE ^ 




FlVlN4 









0.05 



0.1 



0.15 
JafA 



0.2 



0.25 



FIG. 11: Phase diagram as in Fig. [5] including the local 
nearest- and next-neighbor OO interactions with k = 0.2t 
and k' = O.lt (spring constants K2/K1 = 0.25). '*' indicate 
representative points in the CE region for which the correla- 
tion functions are presented here. The dashed arrow presents 
the direction of the transition from the localized CE to the 
itinerant electron A-phase at low temperatures. 



where — (+) refers to i <E A (i 6 B) sublattice, respec- 
tively, and \4>\ < j. Here, \ix) and \iz) stands for local 
basis orbitals: |a;) and \z) at site i, respectively. In the 
rotated basis {\ifJ-}, we assume for the transformed 
operators: (T/) = -1/2 (1/2) for i S A (B) sublattice, 
respectively. At intermediate values of 4> one can reach, 
e.g., \3x 2 -r 2 )/\3y 2 -r 2 ) (\z 2 -x 2 )/\z 2 -y 2 ))-type orbital 



ordering at <f> = 7r/12 (<fr = —tt/12), respectively. The ro- 
tation (|l(jf) does not include complex orbital order which 
was proposed as a possible orbital order in manganites 
at smaller doping The classical energy of the distorted 
lattice Ei (per Mn site) is given byjS 



E l = - {Kx + 2K 2 ) [51 + (6' x ) 2 + 51] + K x {u\ + u 2 ) , 

normalized per Mn ion, where Ki (K%) are the nearest- 
neighbor Mn-0 (Mn-Mn) spring constant, respectively. 
Here the force constants include the factor b 2 , and hence 
have the dimension of an energy. 

Next, we have to add the energy of the breathing mode 
(BM) describing the attraction between oxygen ions and 
the unoccupied Mn +4 site: 

E br = ^/3A (5 X + S' x + S z + Ul + u 2 ) ■ (18) 

Here, all oxygens are equally attracted towards the cen- 
tral manganese ion. The BM induces charge alternation 
on neighboring Mn sites and contributes to the effective 
inter-site charge-charge repulsion (V). The coupling con- 
stant to the breathing distortion is written as /3A, where 
A denotes the coupling to the JT distortion. The term 
describing the JT coupling depends on the character of 
the occupied orbital at a Mn +3 site. The JT interaction 
describes the coupling between the occupied e g orbital 
at a particular Mn 3+ ion and the distortions of the sur- 
rounding oxygen ions, as introduced by MillisS In the 
orbitally rotated state l|16f) this energy has the following 
form, 



1 



^Jt(^) = 77 A 



S z + iti + it 2 - ^ (S x + S x ) 



y3cos(2^) 



- (S x - 5 X ) + u 2 



sin(2</>) 
(19) 



and depends on the rotation angle <fi and on the defor- 
mation of the lattice. 

The optimal values of the ionic displacements for a 
given angle cj> are found by minimization of the total en- 
ergy, E tat (<t>) = Ei + E br + S JT (</•), which for (T7J)-® 
and neglecting higher-order nonlinear terms gives, 
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This yields for the total energy of the distorted lattice, 
A 2 



EtM = — 



{f3 2 (5K 1 +4K 2 



4,K l {K 1 + 2K 2 ) 
+ {Ki + 2K 2 ) [2^sin(20) + cos 2 (20)] + 2K X + K 2 } . 

(21) 

This implies a change of the volume of the lattice with 
respect to the ideal perovskite (Vo — b 3 ) SV/Vo ~ 8 X + 
S' X +5 Z = —3 A/3/ (K\ + 2K 2 ). Hence the volume decreases 
with increasing BM coupling (3, but does not depend on 
the character of the orbital order. 

The relative strength of the effective next- and nearest- 
neighbor JT interactions can be obtained by comparing 
the lattice stiffnesses of the orbital modulations in the 
undoped (~ k) and the quarter-filled (~ k') case, re- 
spectively. For undoped LaMn03 the lattice energy was 
found to be [see Eq. (2.14) in Ref-EU, 



3A 2 



K 1 {K 1 + 2K 2 ) 



{K! + K 2 [l + cos(40)]}. 



(22) 

In the half-doped case the orbital alternation at Mn 3+ 
neighbors is driven by the Mn 4+ ion shift. The energy 
(21) can be decomposed E to t{4>) — Eo(4>) + 8E((j)), where 
5E(<fi) is the gain of energy due to the Mn 4+ displace- 
ment: 



6E(<P) = 



3A 2 



lQ{Ki + 2K 2 ) 



[1 + 008(40)] . (23) 



This cooperative energy does not depend on the breath- 
ing mode coupling (3 directly, though there is an indirect 
dependence via Eq(4>) influencing the optimal value for 
(j). Contrary to the undoped case (|2*2*|l . where the Mn- 
Mn potential (~ K 2 ) is necessary to obtain finite orbital 
stiffness, in the half-doped case I)2.'{[1 the Mn-0 interac- 
tion alone can lead to cooperative JT effect stabilizing the 
alternating orbital order. Assuming A = 6eV, K\ = 200 
eV, and < K 2 jK\ < 1 (see Ref . I57ll62|l one finds for the 
|x)±|js) ((f) = 0) orbital order the CE energy contribution 
per Mn ion 8E ~ 20 - 70 meV. 

In Table |U we present the ratio 8E((f)/Etot{<t>), charac- 
terizing the stability of the CE orbital order, as function 
of the BM coupling strength j3 and the ratio of spring 
constants K 2 jK\. The largest relative gain in energy 
due to the Mn 4+ ions shift (~ \8 X — S' x \b/2) is obtained 
for the (\x) + \z))/(\x) — \z)) orbital order (<j> = 0) which 
for small BM coupling (~ (3) can reach up to 50% of 
the total lattice energy, whereas for larger values of (3, 
although 6E(<f>) is not changed, the total lattice energy 
increases rapidly as ~ f3 2 being dominated by the dis- 
tortions of the MnOe octahedra rather then by their 
displacements. Furthermore, with increasing Mn lattice 
stiffness (~ K 2 /K\) both the 8E(4>) / E to t(4>) ratio and 
the shift of unoccupied Mn 4+ ions [~ A/ (K\ + 2K 2 )\ de- 
crease. 

Comparing the cooperative part [~ cos(40)] of both 
JT energy contributions [ (|22|l and l|23|l ] we end up 



TABLE I: The relative energy gain, 8E(cj>)/ 'E to t(4>), associ- 
ated with the Mn 4+ ion displacement (see Fig. IlUt for the 
orbital orderings given by (f> = and ±7r/12 presented as 
function of j3 and K2/K1. 
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FIG. 12: Correlation functions in the case of a CE ground 
state (as in Fig. |HJ but obtained for the local OO interaction 
(26) with spring constants K2/K1 = 0.25 (k = k/2). Other 
parameters: V — t, k — 0.2t, Jaf = 0.04t, Jh = 15t. 
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FIG. 13: Charge correlations in the case of a FM ground state 
as obtained for the local OO interaction with spring constants 
K 2 /Ki = 0.25 (k = k/2). Other parameters: V = 0.5*, 
k = 0.2t, Jap = 0.02*, Jh = 15*. 



with a simple relation between different OO interaction 
strengths, 



K 
K 



8K0 



(24) 



For K 2 = Ki the ~ i? -3 scaling of elastic interactions^ 
is recovered. In the above estimate we have neglected 
the superexchange interaction which can also contribute 
to the effective k [see e.g. Eq. (2.7) in Ref . E3 - 

The relative strength of both spring constants in the 
CO case can be assessed from the relative length of differ- 
ent bonds in the half-doped Mn02 plane. As measured 
for low-temperature Lao.sCao.sMnOa superstructure 4 ^ 
the length of all four in-plane Mn 4+ -0 bonds is very sim- 
ilar (d 2 = d 3 ~ 1.9151). Neglecting the effect of Mn-O- 
Mn bond bending and assuming d 2 — ^3 (see Fig. HO)) 
one obtains the relation, 



(25) 



by using Eqs. (J2DJ. The above assessment together with 
(|24|l underlines the importance of Mn 3+ -Mn 3+ orbital 
coupling (~ «') in CO manganites. 

The structural data of Radaelli et al~ may be used 
in combination with Eqs. (|20l a-e) to infer the values 
for A and (3. We have not attempted to determine the 
force constants as well, but assume 5 - 7 - K\ = 200 eV . 
Using basically data for the JT-distortion of Mn 3+ -0 
in-plane bonds (di ~ 2.07 A and 1I4 ~ 1.92 A) and the 
Mn 4+ displacement, we obtain for K 2 — 0.25(0.0)^1 
/3 ~ 0.27(0.35) and A ~ 7.4(5.7) eV, respectively. We 
note that the ratio \/K\ ~ 0.037(0.029) is rather close 
to the values determined by Millie 57 for LaMnOa. This 
also implies that the energy gain due to the breathing 
distortion is substantial [see Eq. (|2ip]. 



E. Effect of cooperative Mn 3+ -Mn 3+ JT 
interactions 

We return to the many-body theory for a translation 
invariant system. As the OO interaction between further 
Mn 3+ neighbors is mediated by an Mn 4+ ion we intro- 
duce the following effective three-site interaction, 



^OO — %K 



<ijj'> 



(26) 



where the operator T\y is defined as in Eq. (|12f> . and (ijj') 
denotes three neighboring sites along a or b direction. 
The reasoning here is similar as for the nearest-neighbor 
JT interaction in Eq. I|llfl . i.e., the expectation value 
of the operator Ty in an orbital ordered state depends 
on the orbital angle precisely as in Eq. 112,'jll . namely ~ 
[1 + cob(40)]. The (1 - nj) factor in Eq. ® reflects 
the charge state of Mn 4+ and favors CO in the quarter- 
filled system, while the contribution of the interaction 
(|26|l would vanish in the undoped compound. 

Investigating the model with both nearest-neighbor 
(fTTf) and three-site OO couplings we find a strong 
change in the low-temperature phase diagram (see 
Fig. II If) as compared with the previous ones (Figs. |3 
and UJ. Assuming k' = k/2 (which corresponds to 
K2/K1 = 1/4), the CE-like correlations set in already 
at rather small inter-site repulsion (V ~ 0.5*) while for 
stronger Coulomb interactions (V > 1.5*) the spin corre- 
lations have predominantly AF character. Furthermore, 
the FM region is quite small and can be found only for 
V < * and Jaf ^ 0.1*. When decreasing the ratio K 2 jK\ 
(and keeping the value of k' constant) the FM region 
shrinks at the expense of the CE region. 

Next, we present the two-site correlations for one rep- 
resentative (V — *, Jaf = 0.04*) point in the CE re- 
gion in Fig. ^] From Fig. ^] (a) one sees that the CE 
state can exist for less pronounced CO than in the elas- 
tic strain model considered in Sec. IIII CI where almost 
perfect Mn 3+ /Mn 4+ CO was necessary to stabilize the 
zigzag state (see Fig. E). Moreover, the evolution of CE- 
type orbital order with decreasing temperature is here 
correlated with the onset of CO [compare Fig. 1121 (a)-(b) 
with|S](a)-(b)]. Regarding (5f mS'f, R ni+R) correlations, 
the CE-type spin order develops at lower temperatures 
than the charge and orbital order but at higher temper- 
ature than previously found for the strain induced in- 
teractions. The three-site term (|26|l which is active for 
a Mn 3+ -Mn 4+ -Mn 3+ sequence along the a and b direc- 
tion, respectively, promotes not only orbital alternation 
on neighboring Mn 3+ sites but also strengthens CO. This 
is clearly seen in the charge correlations which are now 
stronger than those e.g. in Fig. [21 although the same 
Coulomb repulsion (V — *) was assumed in both cases. 

Going from the CE to the FM phase (as indicated by an 
arrow in Fig. Ill|) one finds that the CDW is now reduced 
at low temperatures in the FM state fFig. H3|) . while the 
charge modulation is still rather strong at higher temper- 
atures (T/t w 0.5) above the Neel temperature where the 
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FIG. 14: Energy gained by e 9 electron stacked in the c direc- 
tion as function of charge modulation S obtained for V = 0.5t, 
k! = O.lt and different values of k. 



electron kinetic energy is quenched due to the random ti g 
spin orientations. This behavior is a clear manifestation 
of the DE mechanism. We note that the CE-type orbital 
correlations are still present in the FM phase in Fig. ^] 
for V > 0.5£, whereas for smaller values for V x 2 -y 2 or- 
bital correlations dominate, which are characteristic for 
the metallic A-phase. 



F. Charge modulation in the c direction 

Finally, we consider possible mechanisms stabilizing 
the charge stacking of CE-ordered planes^ as observed 
in most charge-ordered manganites at doping x = 1/2. 
Any interaction stabilizing such a stacking must over- 
come the inter-site charge-charge repulsion. There are 
two interactions which contribute to such a repulsion: 
(i) the electron-electron Coulomb force V; and (ii) the 
BM coupling. Both interactions will effectively tend to 
maximize the number of Mn 3+ -Mn 4+ bonds. The BM 
coupling contribution to the inter-site charge-charge re- 
pulsion (Vbm) can be estimated comparing the lattice 
energies with charge alternating and stacked along the 
c direction which leads to the estimate Vbm = A 2 [/3 — 
sin(2cj))] 2 /(2Ki). Assu ming typical values for A = lOi 
and K\ = 50CM (see Ref. l57l) together with our previous 
estimate (3 ~ 0.35 one can obtain e.g. Vbm — 0.012£ for 
the {(\x) ± \z))/y/2} CE orbital order. 

In a multi-layer system the 00 interaction along the c 
direction is given by the nearest-neighbor coupling ~ k, 
which in the {\x}, \z)} basis has the form^ 



H, 



oo 



(27) 



This interaction together with the ~ k' coupling, which 
controls the CE-orbital order in the (a, b) plane, in a 
charge-ordered case yields an approximate orbital model 



oo 



6k' 



E 



T z T z 



(28) 



where T? (Tj x ) are dehned with respect to the {(\x) ± 
\z))/V2} orbital basis. Defining |±) = Qx) ± \z))/y/2 
and assuming that in the (a, b) planes of a bilayer sys- 
tem and |+)j orbitals are occupied on two neigh- 
boring sites in the c direction, the ~ k coupling [see 
Eq. J2HJ] rotates orbitals towards the |— )i<8 |— )j configu- 
ration and leads to binding of two electrons with energy 
Ef, = \J (12k') 2 + Ak 2 — 12k'. For charge disproportion 
with Mn 3 ' 5 * 5 ions alternating in the (a, b) planes and as- 
suming that the k (k') two (three)-site interaction scales 
as ~ (1/2 + 8) a with a — 2 (3), respectively, one finds 
for the energy gain (per Mn ion) due to charge stacking 
along c direction for a 3D lattice: 



1 



SE = 6k'(- + S) 



36K' 2 (i+<5) 2 



1 



2S 2 V. 



(29) 

The last term stands for the effective Coulomb repulsion 
at given charge modulation. With charge modulation 
sufficiently reduced S - 0.2 - 0.3 (see Fig. EJl the OO 
inter-plane coupling can stabilize charge-stacked phase 
(SE > 0). Such stabilization is strengthened (weakened) 
by increasing k (k') OO interaction, respectively. For 
smaller charge disproportion (5 < 0.1) and considerable 
charge on Mn 3 5-5 sites the ~ k in-plane coupling [ne- 
glected in Eq. J2EJ] can also play an important role in the 
inter-plane charge modulation. However, when the CO is 
strongly reduced the orbital ground state is too complex 
to be analyzed analytically. 

A further mechanism stabilizing the observed e g charge 
stacking along the c direction results from the AF inter- 
plane superexchange interactioniSS Yet, the magnetic in- 
teraction (Jaf) results mainly from the tig — t2g su- 
perexchange which is of the order of only a few meVi^ 
Thus, the antiferromagnetism in the c direction observed 
in most of the CE structures appears to be a result 
rather than the origin of the stacking pattern of the 
charge /orbital order. 



IV. SUMMARY 

We have investigated the potential of different electron- 
lattice mechanisms to stabilize the CE spin-orbital or- 
dered phase in charge ordered half-doped manganites 
(filled with one e g electron per two Mn ions). Our study 
is based on a generalized DE model containing the e g 
orbital degrees of freedom, and also including inter-site 
Coulomb repulsion (V) and AF superexchange (Jaf) be- 
tween core spins. The effect of three types of orbital 
interactions derived from the cooperative JT effect have 
been studied: (i) nearest neighbor JT-coupling, which is 
responsible for the alternating orbital order in undoped 
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LaMnC>3, (ii) in addition second neighbor JT- interaction 
with a strength as expected for elastic strain, and (hi) 
a 2nd neighbor JT interaction which accounts for the 
breathing distortion of the Mn 4+ octahedra and is in gen- 
eral stronger than (ii). A particular feature of this inter- 
action, which has been proposed based on the classical 
model considerations in Sec. lIIIDl are the displacements 
of the Mn 4+ ions in the CE-phase. 

We have studied these models employing finite tem- 
perature diagonalization on a v° x y8 cluster with pe- 
riodic boundary conditions, thereby keeping the transla- 
tional invariance. The evolution of the different phases 
with temperature is monitored by the calculation of cor- 
relation functions for charge, orbital and spin, respec- 
tively. Starting out from the generalized FM Kondo lat- 
tice model Q without any orbital interaction (Hoo = 0) 
we find the orbital correlations dominated by the x 2 — y 2 
orbitals which are favored by the kinetic energy42*2i This 
orbital state is characteristic for the metallic A-phase 
in wide-band manganites at half doping. For small AF 
superexchange interactions (Jaf <C t) between ti g core 
spins the FM state is found in the spin sector, while for 
larger values of Jaf the ground state is that of a regular 
2 sublattice antiferromagnet. Even for large Coulomb 
repulsion V we found no evidence for C or CE struc- 
tures for some intermediate values of the parameter Jaf- 
The inclusion of the nearest-neighbor JT- interaction (i), 
however, was found to stabilize the AF-C phase in a wide 
parameter range in the V — Jaf phase diagram. 

To obtain the more complex CE orbital and spin state 
not only the kinetic energy must be reduced by increas- 
ing the Coulomb repulsion V but also a further neigh- 
bor antiferro- orbital interaction must be included in the 
model. This interaction results from the electron-lattice 
coupling and is a generalization of the JT interaction 
which was found to play an important role in the parent 
LaMnC>3 compoundSS In CO compounds at quarter- 
filling the effect of the 00 coupling will be dominated by 
interactions between further distant Mn neighbors. We 
have shown that such interactions can produce the CE 
structure in both spin and orbital sectors. In the case 
of the three-site effective interaction (~ k') the zigzag 
chains are found to be stable not only for modest inter- 
site repulsion (V ~ 0.5i) but even in the limit of V — 
and provided k is small, where the CO is induced en- 
tirely by the three-site orbital coupling (|26|l . This effec- 
tive next-nearest-neighbor 00 coupling was derived tak- 
ing into account the cooperative nature of both Mn 3+ -0 
and Mn 3+ -Mn 4+ bonds length deformation. The free 
relaxation of the latter length leads in presence of CO 
to orbital alternation in planar directions, and the pat- 
tern of directed orbitals typical for the orbital CE struc- 
ture appears. Consequently, with the help of the orbital 
topology; 1 ^ 3 - this gives rise to the CE AF spin ordering 
for realistic values for Jaf- 

An important aspect of the second neighbor JT- 
interaction (iii) is its 3-site structure which involves a 
coupling to the charge at the Mn 4+ ion. This implies 



that the evolution of the CE orbital correlations with 
temperature are linked to the evolution of CO, while the 
magnetic CE correlations develop at lower temperature 
similar to Pr 1 / 2 Cai/2MnO J 36 i 37 . This is different from 
(ii), where first CO develops, then 00 and finally at 
lowest temperature AF spin order. 

As discussed in Sec. IIII Al we did not find the CE spin 
and/or orbital phase considering a pure electronic model, 
i.e., neglecting the electron-lattice coupling. Although, 
the topological effect combined with the DE mechanism 
leads to the preference of the CE over C phase when 
all e g (bridge and corner) orbitals are degenerate^ 3 - the 
situation is more complex when such degeneracy is lifted 
by JT lattice distortions. As described in the Appendix, 
one finds band insulators in both cases while the stability 
of the CE versus the C phase strongly depends on the 
crystal- field splittings (E z , r\E z ~ A). 

An important issue is the degree of charge ordering 
in the CE phase. The CE spin-orbital structures ob- 
tained in our simulations appear in the regime of rela- 
tively strong CO, as infered from our calculated charge 
correlation functions. Our calculations also suggest that 
the second neighbor JT coupling even for small or vanish- 
ing Coulomb repulsion can lead to a pronounced charge 
modulation. These results are consistent with the simple 
physical picture, that DE perpendicular to the FM chains 
must be strongly suppressed to allow AF superexchange 
to trigger the CE structure, which may not be the case for 
weak CO. On the other hand, a recent study by Mahade- 
van, Terakura and SarmaS& based on the band structure 
approach arrived at the conclusion that the charge dif- 
ference between the two Mn species in the CE phase is 
basically negligible, which would imply the absence of a 
significant breathing distortion of Mn 4+ octahedra. Fur- 
ther experimental information on the degree of lattice 
distortion in the CE phase will certainly help to test the 
proposed JT based mechanism for the stability of the CE 
phase. 

The charge modulations determined in our calculations 
may be influenced by numerical limitations which overes- 
timate quantum fluctuations: (i) spins S = 1/2 are used 
to describe t 2g electrons; (ii) the calculations are made 
using a 2D cluster; (iii) the oxygen ions are not present 
explicitly in the cluster and thus only dynamical interac- 
tions are included (k, k'), which are mediated by oxygen 
displacements, while possible static lattice effects^! are 
neglected. Certainly each of these approximations weak- 
ens the charge ordering and the spin-orbital order. It is, 
however, not straightforward to decide whether the re- 
gion of the CE-phase in the phase diagram may grow or 
shrink, if quantum fluctuations are suppressed. 

Finally we note that further neighbor JT interactions 
can also play an important role in the charge stripe for- 
mation at other commensurate carrier concentrations. 
For x < 1/2 (x > 1/2) the dominant coupling arises from 
the nearest- (next nearest-) neighbor 00 interactions, re- 
spectively, with a pronounced electron-hole asymmetry 
present. In the latter case with small inter-site Coulomb 
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repulsion (V) the bi-stripe stated would be preferred 
while larger values of V would favor a "Wigner crystal" 
charge arrangement ^ 
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APPENDIX: ID BAND MODEL 



a, o.o ■ 
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A quite appealing explanation for the stability of 
the CEas compared to the C phase was suggested in 
Refs. Illlbil by considering a ID model, where in the CE 
phase a gap opens as a consequence of a topological sign 
in the hopping matrix-elements. However, the intra-site 
Coulomb repulsion U (leading to some redistribution of 
the e g charge) can destabilize the CE phase. 6 -'—. Here 
we investigate how these arguments are modified when 
orbital degeneracy is lifted due to JT distortions. We 
consider the Hamiltonian: 



H 1B = H t 



E 



z(b) 



E 



T?+\riE <b) 



E «*« (A-l) 

i£bridge 



which contains the kinetic energy H t of an electron hop- 
ping along a FM chain 13 with possible e g level shift at 
the bridge sites (oc rjE z Q,\) and level splitting at corner 
sites (oc ) considered in two different orbital basis 
{!*>, \z)} and {\a) = -^(\x) + |6> = j-{\x) - \z))}. 



We assume that the energies of the bridge 3x 



r 2 orbitals are located between the split energy levels at 
the corners assuming |7j| < 1. Both kinds of level shifts 
can result from the lattice distortions. At the bridge site 
only the orbital 3x 2 — r 2 (3y 2 — r 2 ) is considered while the 
orthogonal one y 2 — z 2 (x 2 — z 2 ), respectively, is decou- 
pled in chain direction. In momentum space this leads 
to a 3 x 3 matrix problem: 



H 1D = 




K 



K 

E z 






-E z 




where Wk = 



(A.2) 

-2tcosk (2tcosk) and W' k = -2\/3t sin k 
{2^/itcosk) are the hopping elements for the CE (C) 
phase, respectively, in the {|x), \z)} orbital basis. 

Calculating the band energies (see Fig.ll5fl one can eas- 
ily determine the character of the stable phase. For the 
degenerate band model (E z — 0) the C phase is metallic 
while in the CE phase a gap opens due to the topological 



FIG. 15: Electron dispersion for the zigzag chain of the CE 
phase (solid lines) and for the linear chain of the C phase 
(dashed lines) calculated for different crystal-field splittings 
in the {\x), \z)} orbital basis at corner sites (E z ). The orbitals 
at the bridge sites are degenerate (assuming 77 = 0). 
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FIG. 16: Phase diagram of the chain model illustrating the 
stability of the CE against the C phase as function of the 
crystal-field splittings defined with respect to the orbital basis 
(a) {|x), jz}} and (b) {-j={\x) ± |z))}, respectively. 



signm^i^ The opening of the gap stabilizes the CE phase. 
When the degeneracy between \x) and \z) states at the 
corner sites (E z ^ 0) is lifted, the C phase becomes insu- 
lating as well and for sufficiently large e g level splittings 
it can have a lower energy than the CE phase when the 
chain is doped with one electron per two Mn sites [see 
Fig. I15f b)]. The region of stability of the CE phase in- 
creases for positive (negative) E z when r\ — > —1(1), re- 
spectively [see Fig- llCf a)]. Generally, the CE phase is sta- 
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bilized with respect to C, when the occupied bridge state 
approaches the lower energy level of the corner states. In 
a real half-doped crystal where in the CE phase Mn— O 
bond lengths shrink (expand) in the c direction [(a, b) 
plane], respectively, one would expect an energetic pref- 
erence, due to the JT effect, for orbitals with stronger 
x 2 — y 2 than 3z 2 — r 2 character. Thus, in the ID model 
the relevant part of the phase diagram (Fig. I16|l is the 
region with E z < and r]E z < 0, i.e., where the CE 
structure dominates. 

To relate our mean-field considerations to the model 
from Sec. IIII Bl the stability of the CE phase versus C 
was also considered with level splitting at the corner sites 
defined with respect to the orbital basis {4=(|ir) ± |z))}. 



Such a splitting, directly related to the nearest-neighbor 
JT interaction 1)11(1 when Eb > 0, is operative in the C 
phase but vanishes by symmetry in the CE orbital phase. 
As shown in Fig. Iltif b) the CE phase is stable for small 
splittings (0 < E b < 2t) and in the E b > 0, rj —> -1 
region. The latter is the relevant parameter range in 
the case of nearest neighbor JT interactions, and hence 
the model suggests the CE phase as stable phase. This, 
however, is in conflict with our numerical study of the 2D 
model, where only the C phase was found to be stable. As 
in 2D the nearest neighbor JT interaction stabilizes the 
C structure, whereas this interaction is frustrated in the 
CE phase. This effect is not included in the ID model. 
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